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ABSTRACT 


Despite their importance for air-sea interaction and microwave remote sensing of the ocean surface, 
intrinsic properties of short gravity-capillary waves are not well established. This is largely due to 
water surface currents and their effects on the direct measurements of wave parameters conducted at 
a fixed point. Frequencies of small scale waves propagating on a surface which itself is in motion, are 
subject to Doppler shifts. Hence, the high frequency tail of the wave spectra obtained from such tem- 
poral observations is smeared. Conversion of this smeared measured-frequency spectra to intrinsic- 
frequency (or wavenumber) spectra requires corrections for the Doppler shifts. Such attempts in the 
past have not been very successful in particular when field data were used. This becomes evident if 
the amplitude modulation of short waves by underlying long waves is considered. Microwave radar 
studies show that the amplitude of a short wave component attains its maximum value near the crests 
and its minimum in the troughs of the long waves. Doppler-shifted wave data yield similar results 
but much larger in modulation magnitude, as expected. In general, Doppler shift corrections reduce 
the modulation magnitude. Overcorrection may result in a negligible modulation or even in a strong 
modulation with the maximum amplitude in the wave troughs. The latter situation is clearly contra- 
dictory to our visual observations as well as the radar results and imply that the advection by currents 
is overestimated. In this study, a differential-advection approach is used in which small scale waves 
are advected by the currents evaluated not at the free surface, but at a depth proportional to their 
wavelengths. Applicability of this approach is verified by the excellent agreement in phase and mag- 
nitude of short-wave modulation between results based on radar and on wave-gauge measurements 
conducted on a lake. 


INTRODUCTION 

Spatial measurements can yield directly the wavenumber spectrum of short gravity-capillary waves 
which is the key quantity in characterization of waves and microwave remote sensing of the ocean 
surface. Spatial measurements are difficult to obtain, so attempts have been made to determine the 
wavenumber spectrum from the frequency spectrum of wave height or slope observed at a fixed 
point in the laboratory (Sinitsyn et al, 1973; Reece, 1978) or in the field (Evans and Shemdin, 1980; 
Ataktiirk, 1984; Stolte, 1984, 1989; Richter and Rosenthal, 1986; Ataktiirk and Katsaros, 1987). 
The approach requires relating the wavenumber and the frequency of wave components via a disper- 
sion relation. In the presence of long waves and surface currents, the dispersion relation for short 
waves is greatly altered (Kitaigorodskii et al, 1975) and the problem is complicated. 

In a Eulerian frame of reference, the observed frequency of a short wave depends on its translational 
speed as it passes by the probe. If a wave component is travelling at its phase speed, c(a))=co/k, then 
the observed frequency is its intrinsic frequency, (o, which is related to the wavenumber, k, via the 
dispersion relation, (o=(g+Tk/g w )k 2 , where g is the gravitational acceleration, g w is the density of 
water, and T is the surface tension. If the waves are riding on a surface which itself is in motion, then 
their observed translational speed is the sum of their phase speed and the component of the surface 
current, U s , in the direction of wave propagation specified by the unit vector k/k. In this case, the 
observed frequency is o=co+k .U s =co+kU s cos<f>. 
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The wave field considered in this study is generated purely by the local winds over the limited fetch 
of a lake, where tide and swell are not present. Then, the surface velocity can be described as, U s (t) = 
U t + Ust + U w , where the three terms on the right hand side are the horizontal orbital velocity of the 
long gravity waves, wave-induced mass transport (Stokes drift), and wind-induced surface drift, 
respectively. Then, the relation between the frequency of encounter, a, and the intrinsic frequency, 
a), becomes; o=co +k . (U^ + Ust + U w ), which is an explicit form of the Doppler frequency shift. The 
dispersion relation appropriate for this case may be described as before with g replaced by a modified 
gravitational acceleration, g', which includes the additional effects caused by the underlying long 
waves (Garrett and Smith, 1976; Phillips, 1981; Longuet-Higgins, 1985). 

In the past, Doppler shift corrections have been made using the currents evaluated at the free surface. 
In the present study it is shown that this approach may lead to unrealistic results. Therefore, an alter- 
native approach in which the short waves are assumed to be adverted by the current evaluated at a 
depth proportional to their wavelengths is tried. The appropriateness of the latter approach is verified 
by experimental results. 

EXPERIMENTAL SET-UP AND INITIAL DATA PROCESSING 

The field data used in this study were collected during the summers of 1986 through 1989 at Sand 
Point, Lake Washington, Seattle. The research facility is operated by the University of Washington, 
Department of Atmospheric Sciences . It consists of a mast located 15m offshore, a hut situated on a 
beached barge and a floating dock between the mast and the land. On the beach, one end of the dock 
is mounted on a rotating base so that the other end can swing between the barge and the mast. In order 
to prevent interference with the measurements, during data acquisition the dock is secured along the 
shore away from the mast. A data link between the instruments on the mast and the data acquisition 
system in the hut is supplied by underwater cables. 

Ideal experimental conditions at this site are achieved during northerly winds when the fetch over the 
lake reaches a maximum of about 7 km. The location offers a unique opportunity to study surface 
waves generated by the local wind on a natural body of water under a variety of environmental condi- 
tions without the complexities that may arise from the presence of swell or tidal currents. The water 
depth, D, by the mast is approximately 4m and increases rapidly further offshore. We estimate the 
degree of reflection of the gravity waves from the inclined beach to be 5% or less (Miche, 1951). 
When kD>n/2, the waves are said to be in deep water (Phillips, 1977, p. 37). This corresponds to 
A^2rc/k<16 m, where X is the wavelength. Since this condition was met in all cases (typically X= 5 m 
for the dominant gravity waves), the location can be considered as representative of deep water. 

During the experiments, the following variables (Table 1) were directly measured or observed using 
the indicated instrumentation. Radar measurements were conducted by the University of Kansas 
group and are described by Salam et al (1991) . Details of our work including various corrections to 
measured signals and data processing techniques are provided by Atakttirk (1984, 1991) and 
Atakturk and Katsaros (1987, 1989). 


Table 1: List of directly observed variables and instrumentation used. 
Radar measurements were conducted by the University of Kansas. 


Wave height 

Resistance wire gauge 

Wind components 

K-Gill anemometer 

Air temperature (dry & wet bulb) 

Thermocouples 

Water temperature 

Thermistors 

Visible surface features 

Video camera i 

Electromagnetic backscatter 

Microwave radar 
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Initial data processing involved the following calculations. Surface fluxes of momentum, sensible 
heat and latent heat were determined by the eddy correlation method for every 1 7-minute segment of 
data. This provided us With the surface friction velocity, w*, atmospheric stratification parameter, 
z/L, and surface roughness length parameter, Zo . Further, a digital linear-phase, band-pass filter was 
applied to the wave height time series, £(t), to separate its high-frequency (3-11 Hz) part, £'(t). A 
Fast Fourier Transform algorithm and a sliding window of 0.5 second long were applied to £'(t) to 
create a time series of the measured frequency spectra, E(a) with a sampling rate of 8 Hz. 

Surface currents were not directly measured but estimated from u* or £(t) via a linear wave theory. 
For a monochromatic wave train, these relationships can be described by 

t,(t ) = A cos(kx — cot) ; U £ = £cd e kz ; U St = (ok A 2 e 2kz ; U w = 0.60 «* (1) 

In the case of field data, contributions from wave components with frequencies up to three times the 
peak frequency were included. In previous studies, Ust has been ignored and U w has been approxi- 
mated by 3% of Ujo, wind speed adjusted to 10m height and neutral atmospheric stratification. 

Doppler frequency shifts were corrected in the frequency domain according to the transformation 

E(o) do = E(( 0 ) da = E(k) dk (2) 

by using two different methods to estimate the effective surface currents. In the first method, which 
has been used by all investigators sited at the beginning of this section, surface currents were calcu- 
lated as above and the effective current was taken as that at the free surface, z=0. The second ap- 
proach was different in the sense that the effective wind drift was taken as that at a depth of 
z=-0.044A (Plant and Wright, 1977, 1980) 

U w (z= - 0.044X) = UJ.Z = 0) - (fa) 1/2 £ (3) 


where Zow is roughness length parameter for water, g a is the density of air, and x is the von Karmdn 
constant (also see, Stewart and Joy, 1974; Donelan, 1978; Donelan et al, 1985; Smith, 1990). It 
should be noted that in determining this particular depth, Plant and Wright (1977, 1980) approxi- 
mated Zow withzo, surface roughness length parameter. The same approximation was also used here. 

Evaluating Equations 1 and 3 requires knowledge of u* and Zo- Although these variables were avail- 
able to us through the atmospheric turbulence measurements, they can not be determined with ade- 
quate accuracy over short periods of time due to sampling variability. To overcome this difficulty, a 
method was developed to obtain an estimate of Zo from the measured wave height spectra such that 

^ = 7.2xl0- 3 U, Qj^) 2 ' 1 (4) 

where 'C rms is the root-mean-square wave height, and Om is the measured equilibrium range parame- 
ter (mean spectral amplitude in the range 1.5<(o/to p <3.5 with (tip being the peak frequency; for de- 
tails see, Ataktiirk [1991]). Success of this approach is illustrated in Figures 1 and 2. 

The tool selected for the final evaluation of the experimental results from this study is the modulation 
transfer function (MTF). Following the definitions by Plant ( 1989), modulation of a variable x(t) by 
the underlying long waves, and the coherence function, y 2 , for this modulation were calculated from 


m(f) 


_ £*(/) . 
kx m \ 2 ’ 


72(f) = 


m tw 

ix(/)i 2 \m 2 


(5) 


where X(f) and q(f) are the Fourier transforms of x(t) and g(t), X is the mean of x(t), and * denotes the 
complex conjugate. Since m(f) is a complex quantity, its full description requires knowledge of both 
magnitude and phase. By convention, phase is zero at the long wave crest, and positive ahead of it. 
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EXPERIMENTAL RESULTS 


Radars operating in C (6 GHz), X ( 1 1 GHz) and Ka (34 GHz) bands provided information on MTF for 
Bragg wavelengths about 3 cm or shorter (for results on C and X bands, see Salametal, 1991). On the 
other hand, quality of the wave data for wavelengths shorter than about 5 cm were questionable. 
Therefore, in evaluation of the results we considered other radar studies as well, specifically Wright 
etal (1980) and Plant etal (1983) which provided information in the wavelength ranges of both our 
radar and wave data . For the purpose of this study, the following summary of the radar results is 
sufficient. In the wind speed range of 3 to 7 m/s and for dominant waves with a frequency of 0.5 Hz, 
the magnitude of MTF varies approximately between 5 and 10 decreasing with increasing wind 
speed. Dependence on wind speed is weaker for longer Bragg waves. Maximum modulation occurs 
somewhere between the crest and 30 degrees ahead of it. Coherence between the radar return and 
wave height is strong (y 2 =0.6 ) near the frequency of the dominant waves but rapidly diminishes 
away from the spectral peak. There is little difference between C and X bands. It should be kept in 
mind that MTF obtained from radar measurements has three components: hydrodynamic, tilt and 
range modulation, or m(f) -m' - it+ r, where i = (-1) 1/2 . For our conditions where the radars were 
mounted about 5 m above the surface with incident angle varying between 30 and 60° , ris negligible 
and ill “4 (Plant, 1990). Then, the magnitude of the hydrodynamic modulation, m', which is the 
quantity provided by analysis using the wave data is expected to vary between 3 and 9. 

Results from wave data were obtained using the time series of wave height and short wave spectra in 
Equation 5 . In Figure 3 , the rows correspond to the cases of measured spectra, intrinsic spectra with 
Doppler correction according to Equation 1 (Method A), and intrinsic spectra with Doppler correc- 
tion according to Equations 1 and 3 (Method B). The first column is the maximum coherence squared 
(at the frequency of dominant waves) for each run, The dashed line indicates the level, y 2 =0. 1 below 
which we assumed that the two signals were incoherent. Such points were excluded from MTF anal- 
ysis. The second column is the phase of MTF for all data points with y 2 2: 0. 1 . Mean magnitude of 
MTF for these points are shown in the third column. The vertical bars indicate ± 1 standard deviation. 
Features to be noted in this figure can be summarized as follows. 

MTF magnitude determined from measured spectra is at first fairly uniform across the wavelength 
range, then it slightly increases for X < 15 cm. This is in accordance with the concept that Doppler 
frequency shift poses more serious problems at higher frequencies. Rapid decrease for X<5 cm is 
suspected to be due to noise in the data. The observed increase with increasing wind speed is contrary 
to radar results, and is due to the Doppler frequency shift becoming stronger with the growth of long 
waves. As indicated by the MTF phase, maximum modulation occurs at the wave crests. 

When Doppler shift corrections were made using the currents evaluated at the free surface (Method 
A), resulting MTF had no resemblance to those of the radars except at the longer wavelengths, where 
Doppler shift effects are not significant. In the midrange, both the coherence and modulation dimin- 
ished. Although they reach generally expected values for X < 15 cm, in this region the MTF phase 
becomes ±180 implying that short waves attain their maximum amplitudes in the long wave troughs. 
This is completely against the radar results as well as visual observations, and can be explained by 
overcorrection for Doppler shifts. 

The results based on the approach that short waves are advected by the current at a depth proportional 
to their wavelengths (Method B), are similar to those of radar studies except for X<5 cm where further 
considerations are needed. Dependence of MTF magnitude on wind speed becomes significant at 
shorter wavelengths where the MTF magnitude decreases with wind. Observed MTF phase is about 
zero indicating that short waves attain their maximum amplitudes near the crests of long waves. 
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CONCLUSIONS AND DISCUSSIONS 


Generally it has been assumed that the short gravity-capillary waves are advected by the currents 
evaluated at the free surface. In this study it is showh that the effective advection velocity is rather the 
current evaluated at a depth which is proportional to the lengths of the short waves considered. This 
conclusion is based on the comparisons of the hydrodynamic modulation of short wave amplitudes 
determined through wave and radar measurements conducted on a lake. Even when this differential 
advection approach is used, some conflicts were observed between the radar and wave results for 
wavelengths about 5 cm and shorter. This is possibly due to an assumption made in Doppler shift 
corrections that the wind, surface currents and waves of all scales propagate in the same direction. 
For our experimental site, this assumption is quite appropriate except for small scale waves. Current- 
ly, we are experimenting with the design of an array of wire staffs to obtain the directional wavenum- 
ber spectrum at these wavelengths. 
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FIGURE LEGENDS 


Figure 1: Comparison of the surface roughness lengths determined from: (o) wave measurements 
on Lake Washington; (-) a model particularly for low winds (Lui, Katsaros and Businger, 1979); 
and, (*) atmospheric turbulent measurements during HEXMAX (DeCosmo, 1991), Dashed line is a 
fit to HEXMAX results. 


Figure 2: Comparison of neutral drag coefficients determined from: (o) wave and (+) atmospheric 
turbulent measurements on Lake Washington; and, (*) atmospheric turbulent measurements during 
HEXMAX (DeCosmo, 1991). Dashed line is a fit to HEXMAX results. 


Figure 3: Magnitude and phase of MTF and coherence squared obtained from 
analysis of the wave data against the short-wave length. The rows (from top to 
bottom) correspond to the cases of measured spectra, intrinsic spectra with 
Doppler correction according to Equation 1, and intrinsic spectra with Doppler 
correction according to Equations 1 and 3. Method-B, bottom row, is preferred. 
Symbols in figures indicate different ranges of wind speeds in m/s such that (3-4: 
...+...), (4-5: — x — ) and (5-7: — o— ). The dashed line in the plots of coherence 
squared indicates a value of 0.1. 
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